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Abstract 



We propose a new scenario of non-thermal production of neutralino cold 
dark matter, in which the overproduction problem of lightest supersymmet- 
ric particles (LSPs) in the standard thermal history is naturally solved. The 
mechanism requires a heavy modulus field which decays mainly to ordinary 
particles releasing large entropy to dilute gravitinos produced just after in- 
flation and thermal relics of LSPs. Significant amount of gravitinos are also 
pair-produced at the decay, which subsequently decay into the neutralinos. 
We identify the regions of the parameter space in which the requisite abun- 
dance of the neutralino dark matter is obtained without spoiling the big-bang 
nucleosynthesis by injection of hadronic showers from gravitino decay. The 
neutralino abundance obtained in this mechanism is insensitive to the details 
of the superparticle mass spectrum, unlike the standard thermal abundance. 
We also briefly mention the testability of the scenario in future experiments. 
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I. INTRODUCTION 



Supersymmetry (SUSY) has been a promising candidate for new physics beyond the 
standard model, not only because it is an attractive solution to the naturalness problem 
associated with the electroweak scale, but also because it provides a natural candidate for 
cold dark matter. Indeed, the relic abundance of the lightest supersymmetric particle (LSP), 
which is stable by virtue of the R parity conservation, may naturally fall into the critical 
density with plausible values of its mass and couplings under the standard thermal history 
of the early universe. This concordance was regarded as one of the triumph of particle 
cosmology in the last century. 

Recent progress on the determination of the the cold dark matter (CDM) abundance, 
however, has cast a shadow on the above simplistic view of the standard SUSY cosmology. 
In particular, the value indicated by the WMAP measurement pushes the preferred regions 
of the SUSY parameter to only tiny and special corners of the parameter space. In fact the 
WMAP data imphes the Dark Matter abundance [1] 

ncDuh' = 0.1126l0:°i^i (1) 

at 95 % C.L., where JIcdm represents the current value of the density parameter of the CDM 
component, and h is the Hubble parameter in unit of 100 km/s/Mpc. This is much smaller 
than unity. On the other hand, assuming the standard thermal history of the universe, 
the thermal relic abundance of the neutralinos can be computed in a standard manner. 
Importantly it is known that in generic regions of the SUSY parameter space, the relic 
abundance computed in this way tends to exceed the value the WMAP suggests [2]. The 
right amount of the neutralino abundance is obtained only in tiny and special regions of 
the parameter space where some enhancement mechanisms of neutralino annihilation take 
place. They include 1) light sfermions, 2) coannihilation with e.g. chargino, stau, or stop, 
3) resonance enhancement, and 4) annihilation into W pair. 

Furthermore this tension we now have may become real conflict if the LHC experiments, 
for instance, will find SUSY particles outside the regions where the calculated relic abun- 
dance is consistent with that implied by the WMAP measurement. 

At this point, one should recognize that the argument of the thermal relic abundance 
strongly rehes on the standard thermal history of the universe. The fact that the enhanced 
annihilation of neutralinos is required to meet observation may imply that the universe might 
have evolved in a non-standard manner. Alternatively, we may find a natural solution 
to the dark matter problem making use of a non-standard thermal history without any 
enhancement of the annihilation processes, which is the central topic of the present paper. 

Tension in the SUSY cosmology with the standard thermal history increases even more 
when one turns to the recent reanalysis of the gravitino problem which includes the effects of 
hadronic showers on the abundance of the light elements produced by big-bang nucleosyn- 
thesis (BBN) [3]. It is well-known that decays of long-lived particles such as gravitinos would 
spoil the success of the big-bang nucleosynthesis because of the production and destruction 
of the light elements by emitted high-energy photons or hadrons [4-8] . If the gravitino mass 
is just around the weak scale, very severe upper bounds come from the overproduction of the 
^Li abundance and the D abundance. We should also stress that inclusion of the hadronic 
shower also put a rather stringent bound for heavier gravitino whose lifetime is even shorter 
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than 1 sec This is because the hadronic shower produces pions and miclcons which interact 
with ambient nucleons in the thermal bath, changing the neutron to proton ratio and hence 
the "^He abundance [5,9]. Note that the argument is based on the standard thermal history 
of the universe in which there is no dilution or extra production of gravitinos after reheating 
epoch of primordial inflation [10,11]. 

It is, however, quite plausible that the universe did not follow the standard thermal 
history described above. This is particularly the case when there are some moduli fields 
whose masses are not far away from the electroweak scale. The coherent oscillation of the 
moduli fields is likely to dominate the energy density of the universe at some epoch in the 
early universe, followed by rather late decays [12]. If the mass is in the electroweak scale, the 
lifetime of the moduli fields is much longer than 1 sec. meaning that the entropy production 
at the moduli decays takes place with reheat temperature lower than 1 MeV. This would 
upset the success of the BBN. The disaster can be avoided if the moduli masses are somewhat 
larger than the electroweak scale [14-17]. 

In a previous paper [13], we have shown that the gravitino abundance can be regulated to 
an observationally acceptable level if we take entropy production due to moduli decay into 
account. In fact we have shown that it is possible to satisfy the severe constraints imposed 
by its hadronic decay [3] even when the reheat temperature after inflation was as high as 
~ lO^^GeV and the new production mechanism of gravitinos from moduli decay is taken 
into account [13]. 

In this paper, we would like to explore this line of argument and propose a new scenario 
in which the neutralino dark matter is produced in a non-thermal way. In our scenario, the 
neutralinos of the thermal origin are all diluted by the entropy production at the moduli 
decay. The dark matter neutralinos are produced, however, by the decay of gravitinos, 
following the moduli decay to gravitinos. 

It is very non-trivial whether this mechanism survives severe constraints from BBN on 
the gravitino decay which produces hadronic and/or electromagnetic showers. We will show 
that in a range of the mass parameters, the gravitino decay survives the constraint from 
the BBN, whereas the abundance of the neutralino LSPs is in accord with the cosmological 
observations. Thus our model offers a consistent scenario of cosmic evolution in the new 
era of precision cosmology where the cold dark matter abundance is much smaller than 
previously expected and the gravitino problem is much more aggravated. We will also see 
that the scenario requires small hierarchy among moduli, gravitino and soft SUSY breaking 
masses. We will argue how this small hierarchy is naturally realized. 

The organization of the paper is as follows. In the next section we outline the basic 
features of our scenario. In §111 we describe consequences of moduli decay, namely, amount 
of entropy and gravitinos from them. Neutralino abundance from decaying gravitino is 
studied in §IV. Then we compare our yield with the constraints imposed by BBN in §V. 
Finally §VI is devoted to conclusion and discussion. 

II. A SUCCESSFUL SCENARIO WITH SMALL HIERARCHY: OUR PROPOSAL 

In superstring theory and supergravity, there often exist scalar fields which have (almost) 
fiat potential. Examples include moduli fields and the dilaton field in classical string vacua. 
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whose vacuum expectation values (VEVs) determine the size and shape of compactified ex- 
tra dimensions and the string couphng constant, respectively. They also include the field 
which is responsible for supersymmetry breaking in a typical gravity mediated supersym- 
metry breaking scenario (the Polonyi field). They typically have Planck-scale suppressed 
interaction, and the VEVs of the order of the Planck scale. We will generically call these 
fields the moduli fields. 

It is expected that these moduli fields acquire non-fiat potential and thus masses in 
some manner. Their masses depend on the mechanism how the moduli are stabilized. In a 
conventional view [12], the moduli fields acquire superpotential from some effect associated 
with supersymmetry breaking. In this case one naively expects that their masses are of the 
order of the gravitino mass, a typical mass scale of supersymmetry breaking in supergravity. 

Recent development of fiux compactification reveals the possibility that moduli yield 
superpotential when higher-form antisymmetric tensor fields acquire VEVs in compact di- 
rections [18]. The fiux compactification fixes some of the moduli fields, but not all. In 
particular, in the type II superstring the overall modulus whose VEV determines the size of 
the compactified manifold is not stabilized in this way. These remaining moduli may acquire 
potential by other non-perturbative effects, and we may expect that their masses are related 
to the gravitino mass scale. [19,20] 

Thus the naive expectation of the gravity mediation of supersymmetry breaking is that 
the masses of the moduli fields and sparticles in the MSSM sector are similar, and is of 
the order of the gravitino mass. Since the sparticle masses are expected to be in (sub-)TeV 
range, the naive guess leads that the gravitino mass as well as the moduli masses are also 
around this range. In fact, the cosmological moduli problem and the gravitino problem are 
consequences of this naive expectation. 

A closer look at these masses, however, reveals possible deviation of the mass scales in 
one or two orders of magnitude. In the original gravity mediation model by Polonyi, the 
SUSY mass scale is given by hand, and the mass of the Polonyi field is indeed comparable 
to the gravitino mass. In more realistic models, the potential of the moduli fields may 
be generated dynamically. A typical example of this comes from gauge dynamics such as 
gaugino condensation [21]. In this scheme, the hierarchically small SUSY breaking scale is 
achieved by the dimensional transmutation with some numerical coefficient of order 10-100. 
At the same time, the moduli mass which is determined by the second derivative of the 
scalar potential can be larger by the same amount than the gravitino mass which is naively 
expected. In fact, this is the case for the well-known racetrack scenario [22] where multiple 
gaugino condensates take place to obtain a non-trivial minimum of the scalar potential. 
Similar observation has been made in other scenarios, including non-perturbative correction 
to the Kahler potential of the dilaton field [23]. For explicit demonstration, see [24]. The 
moduli masses can also be heavy in some scenarios of fiux compactification [20]. 

This consideration sets a small hierarchy between the moduli masses and the gravitino 
mass. Another small hierarchy between the gravitino mass and the sparticle masses in 
the MSSM can be realized when the field responsible for the SUSY breaking is somewhat 
sequestered from the MSSM sector [25-27]. It is interesting to note that the coupling of the 
moduli fields to the MSSM sector does not violate this property of sequestering, because the 
moduli fields are heavy so that their auxiliary fields have SUSY breaking VEVs suppressed 
by the moduh masses relative to the gravitino mass. Then the softly SUSY breaking masses 
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for sparticles are a combination of the contribution from the moduh SUSY breaking and 
that from the anomaly mediation. [26,28] 

The small mass hierarchy among the moduli, gravitino and sparticle masses given above 
yields a successful scenario without the difficulties of the gravitino as well as the neutrahno 
abundance. Here we would like to sketch our proposal. Prom now on, we assume, for 
simplicity, that there is only one moduhis field which is relevant in our discussion. 

The heavy modulus field is expected to obey damped coherent oscillation with the initial 
amplitude of order the Planck scale. The modulus oscillation dominates the energy density 
of the universe until it decays to reheat the universe. In the present case, since the modulus 
is rather heavy, its lifetime can be much shorter than 1 sec. The entropy production at 
the reheating dilutes the primordial abundances of the gravitinos. Non-negligible amount of 
gravitinos are, however, produced by the modulus decay, on which the constraint from BBN 
will be imposed. 

The sparticles in the MSSM sector are also diluted by the modulus decay, if the decay 
products do not contain the sparticles and the thermal bath after the reheating is not hot 
enough to regenerate the sparticles by scattering processes. We will check that the former 
is indeed the case when the modulus coupling to the MSSM sector is of certain type, and 
the latter will be guaranteed if the reheating temperature is low enough. In this case, the 
neutralinos which remain until today come from the gravitino decay. The main assertion of 
this paper is that in some regions of the parameter space the neutralinos can be abundant 
to constitute the dark matter of the universe, whereas the gravitino abundance satisfies the 
constraint from the BBN. 

The rest of the paper is devoted to describe the details of our proposal given above. 



III. MODULI DECAY: ENTROPY PRODUCTION AND GRAVITINO 

PRODUCTION 

The modulus field is typically displaced from the true minimum during the inflation, 
which starts damped coherent oscillation when the expansion rate of the universe decreases 
to the modulus mass. Unless the initial amplitude of the modulus field (pi is very small 
compared to the Planck scale, the energy density of the modulus field in the form of the 
coherent oscillation will soon dominate the universe. The modulus eventually decays. The 
decay depends on how it couples to other particles [29]. Here we assume that it decays mainly 
to the ordinary particles in the MSSM and thus the universe becomes radiation dominated 
again. This is especially the case when the modulus coupling to the gauge bosons and 
gauginos through gauge kinetic function gives the largest contribution to the decay, where 
the modulus decays dominantly to a gauge boson pair while the decay to a gaugino pair 
suffers chirahty suppression. We write the total decay width of the modulus field 

r„ = iv-^, (2) 

where is the mass of the modulus field, Mq = 2.4 x 10^^ GeV represents the reduced 
Planck scale, and N parameterizes our ignorance of the interaction of the modulus. It is 
natural to expect that ~ 10~^ — 1. Given the above expression, the reheat temperature 
after the modulus decay Tr is calculated to be 
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where g'* counts the effective relativistic degrees of freedom with — 228.75 for the MSSM 
particles. 

Entropy production is also computed in a straightforward way. The entropy increase 
factor is found to be [13] 



A^JA^ A (4) 
V97Vm^ \MgJ ^ ' 

Thus, as far as the modulus mass is much lower than the Planck scale and the initial 
amplitude of the modulus oscillation 0j is not extremely small, the entropy increase factor 
is much larger than unity, which we assume to be the case. The gravitino abundance as well 
as the sparticle abundances before the reheating will be diluted. 

Let us now discuss the production of the gravitinos by the modulus decay. It has been 
pointed out in [30] that the modulus field generically decays to a pair of the gravitinos via 
the following terms in the supergravity Lagrangian: 

Ant = -e^'/^H/^M^^'V^. + H.c. (5) 

where ip^ represents the gravitino field, K is the Kahler potential, and W is the superpo- 
tential. In the above expression, we have the reduced Planck scale Mq to be unity. Since K 
and W have the moduli dependence, the above term in the Lagrangian gives the coupling 
of the modulus to the gravitino pair. Explicitly this is given by 

(e^/2iy) ^ = ^e^'I'^W + e^'I'^W^, (6) 

where the suffix (p here stands for the derivative with respect to 0. Now the auxiliary 
component of the modulus field is written 

As we assume that the modulus field is not a dominant source of SUSY breaking, we can 
approximately equate W<^ = —K^W. Using this, (6) can be rewritten as follows: 

{e^'^w)^ = -\K,e^''W = -\K,m,/,. (8) 
Thus the partial decay width to the gravitino pair becomes of the form 

r(,.2«^c!j™,^cg(!^)\ (9) 

where we have recovered the reduced Planck scale Mq explicitly. In this expression, we have 
introduced a parameter C ~ l/327r|X^/MGp, where a possible numerical factor as well as 
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the phase space factor is ignored. We expect that C ~ 10~^ — 10~^ would be natural, but a 
smaller C may also be possible. 

It has recently been shown that a more effective decay channel to the gravitino and the 
modulino, the fermionic component of the modulus supermultiplet, is kinematically allowed 
if a certain coupling of the modulus to the field responsible to the SUSY breaking is large 
enough [13]. Whether this channel is allowed or not is, however, quite model dependent. 
Here we conservatively consider the former decay mode to the gravitino pair, which will also 
provide a preferred result as we will see shortly. 

Gravitinos are also produced by scattering processes in the thermal bath after the re- 
heating epoch. The production rate is proportional to the reheat temperature and the ratio 
of the number density of thermally produced gravitinos to the entropy density after modulus 
decay reads [31,32,3], 

With the reheat temperature much below 1 TeV which we are interested in, the gravitino 
production by the scattering processes is negligibly small. 

As was shown in [13], the gravitino problem will be solved in this scenario. The primordial 
abundance of the gravitino generated before the reheating due to the modulus decay will be 
diluted to a negligible level. Although the gravitinos are produced by the modulus decay, 
the abundance was shown to be small enough to survive the constraints from the big-bang 
nucleosynthesis in an appropriate region of the parameter space. Here, we consider the decay 
mode of the modulus to the gravitino pair. Then the allowed region in the parameter space 
will become larger than that of Ref. [13]. 

More explicitly, let us estimate the gravitino abundance produced by the modulus decay. 
Here we approximate the decay of the modulus field as an instantaneous decay. The energy 
density of the modulus field is equal to that of the radiation at the reheating, i. e. at the 
modulus decay: 

PR = 1114,714, (11) 

where n,^ is the number density of the modulus. The gravitino number density at the 
modulus decay 113/2 can also be easily estimated in the approximation 

713/2 = 25n^ (12) 

where B represents the branching ratio of the modulus decay into the gravitino pair 

^.B,,^2«4(!^)\ (13, 

The gravitino abundance in the form of yield F3/2, the ratio of the number density relative 
to the entropy density, is given by 

_3 n3/2rii ?> Tr 1 o^,in-io^v-i/2/i7. rn3/2 \V y^/^ 
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IV. NEUTRALINO ABUNDANCE 



In this section, we would like to discuss the rehc abundance of the neutralinos. Here 
we assume that the R-parity is conserved and the hghtest neutralino becomes the lightest 
superparticle, and thus it is a candidate for dark matter of the universe. 

The entropy production at the modulus decay completely dilutes the primordial abun- 
dance of the neutralinos and other superparticles. Thus in our scenario the neutralinos 
which remain today should be produced at and after the reheating due to the modulus 
decay. There are three possible sources of the neutrahnos: 

• Neutralinos produced by gravitino decay 

• Neutralinos produced by modulus decay 

• Neutralinos produced in the thermal bath during and after the reheating 

We now show that the first contribution can dominate over the other two. 

As we discussed earlier, the universe is reheated at the modulus decay. If the reheat 
temperature were too high, the neutralinos would again be in the thermal equilibrium. In this 
case the standard computation of the neutralino relic abundance would apply, resulting in 
the potential problem of the too large relic abundance unless some enhancement mechanism 
of their annihilation is at work. This would be the case if the reheat temperature exceeded 
~ m^/25, where is the mass of the lightest neutralino. On the other hand, if the 
reheat temperature is fairly low, the regeneration of the neutralinos in the thermal bath is 
neghgibly small. Therefore we will consider the case where Tr < m^/30 in the following. 
Note that this condition is not the constraint, but the requisition in the current interest. 

Let us next discuss the neutralinos produced by the gravitino decay. The point here is 
that the gravitino decays eventually to the neutralino LSP under the R-parity conservation, 
and the number of the neutralinos is essentially the same as that of the gravitinos. Thus 
the yield of the neutralino LSPs is evaluated as 

LS'P 3/2 y^^^j VlO^TeVy VlO^TeVy ^ ' 

This is a constant of time, if there is no other process which will change the number of 
the neutralinos nor the total entropy. We would like to emphasize that the neutralino 
abundance obtained in this mechanism is insensitive to the details of the superparticle mass 
spectrum, unlike the standard thermal abundance because Ilsp is directly related with only 
1^3/2- Annihilation of the neutralinos, if it took place, would change it. We will come back 
to this issue later on. On the other hand, using the density parameter of LSP, Ilsp is also 
represented by 

where we adopted the present photon temperature Tq = 2.725 ± 0.002 K [33]. 

Neutralinos may also be produced by the modulus decay. Let us again consider the case 
where the modulus couples to the bilinear of the supersymmetric form of the gauge potential. 
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In this case, the decay to the gaugino pair is known to receive chirahty suppression, and 
thus the branching ratio is roughly estimated as 




Br(0 ^ 2g) - ^ (17) 



where g symbohcally denotes a gaugino. This should be compared with the branching ratio 
to the gravitino pair as 



Br(0 ^ 2g) N ( m-g 



Br(0 2-0^) C \m3/2^ 



(18) 



Note that we are considering the case where {mg/ms/2)'^ ~ 10~^ — 10~^. Thus the contribu- 
tion to the neutralino abundance from the direct modulus decay is sub-dominant, as far as 
the ratio N/C does not exceed ~ (9(10'^). 

The addition of the coupling of the modulus to two chiral multiplets {e.g. quarks and 
their superpartners) does not change the conclusion given above, because the decay to a 
squark pair, for instance, is suppressed by the fourth power of the mass ratio rriq/ 1113/2, 
unlike the second power for the case of the gaugino pair [29]. Thus it is negligible, provided 
that the coupling of the modulus to the gauginos is not small. In the following, we shall not 
consider the direct modulus decay to superparticles in the MSSM. 

In the following we assume that the pair annihilation of the neutralinos after the pro- 
duction by the gravitino decay does not take place. Whether the annihilation is effective 
or not can be seen by comparing the annihilation rate to the Hubble expansion rate of the 
universe when the neutralinos are produced by the gravitino decay. Although the annihila- 
tion cross section of the neutralino iTann is a function of the superparticle mass spectrum, it 
is generically bounded as 

^ Oil ( a\ 

with the coupling cti ~ 10 ^. Then, the ratio of the annihilation rate Tann — ^LSP'^annj 

to 

the Hubble expansion rate i7 ~ T^/ Mq is bounded as 

Therefore, we find Fann/ < 1 at the temperature of gravitino decay. 



V. RESULTS 



We are ready to present our numerical results now. We first discuss constraints we impose 
on the parameter space. Constraints on the reheat temperature after moduli decay are that 
it should be 1) higher than the lower bound required by the successful BEN, Tr > 1.2 
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MeV [16], ^ but 2) lower than the temperature above which the neutralinos are thermalized. 
We set the latter condition as < mLSp/30. 

We show the constraint on gravitino decay from the BBN in Fig. 1. The upper bound 
on the yield of the gravitino 113/2/8 is drawn as a function of the gravitino mass. Effects of 
hadronic shower are taken into account [3] . Here we take the hadronic branching ratio to be 
0.8 because such a heavy gravitino can decay into all of the gauge boson and the gaugino 
pairs and the chiral fermion and its scalar partner pairs. ^ For gravitino mass lighter than 
10 TeV, the constraint comes from the overproduction of ^Li and D abundances. On the 
other hand, the case with heavier mass is constrained by the ^He overproduction. This is 
because the inter-converting reaction between the neutron and the proton becomes effective, 
which is induced by the pions and nuclcons emitted by the hadronic decays of the gravitinos. 
Here we adopted the observational value of the ^He mass fraction Yp{IT) = 0.242± 0.002 
(stat) (± 0.005 (syst)) by Izotov and Thuan (2004) [34]. ^ In the same figure, the dotted 
strip indicates the region where the LSP abundance from the gravitino decay corresponds to 
the desired value Qlsp^^ = 0.0945 - 0.1287 (95% C.L.) in Eq. (16). Here we have assumed 
the LSP mass, mLSP, to be 100 GeV. One can easily scale the strip for other values of ttilsp 
by using Eq. (16). 

By equating Eq. (15) with Eq. (16), we see that the allowed value of is proportional 
to C^/^ for a fixed value of JlLSp/i^- In Fig. 2, we plot the allowed regions for suitable 
observational value of Qi^sph^ = 0.0945 - 0.1287 (95% C.L.) in the C-m^ plane. We plotted 
the strips for 77^3/2 = 50, 10^, 5 x 10^, and 10^ TeV, respectively. We considered only two 
conditions that Tr < mLSp/30 and > 2m3/2- Here we took = 1, and ttilsp = 100 GeV. 
As a result we find that C should lie approximately in the range of 10"*^ ^ C ^ 10~^. Note 
that the theoretically natural value of C ~ 10^^ — 10^^ is completely included in this range. 

We can now identify regions of the parameter space where the LSP neutralino from the 



Here we adopted a milder lower bound on the reheat temperature when we assume that the 
modulus decays only into photons. On the other hand, if we consider the case that the decaying 
modulus also emits hadrons, the lower bound becomes higher Tr > 4-5 MeV [16]. However, this 
difference does not change our results at all for the current purpose. 

^According to the treatments in Ref . [3] , we could not include all of the hadronic decay modes into 
the BBN computation at one time. Therefore, here we consider only the most dominant mode into 
the hadrons, i.e., the gluon-gluino pairs. Then the branching ratio into the gluon-gluino pairs is 
approximately 0.5. Because we adopted this value for the BBN computation, the obtained bound 
on 2/3/2 = ^^'3/2/'^ in Fig- 1 might be milder by ^ 30% approximately. Note, however, that this 
difference of the hadronic branching ratio does not change the constraints on the model parameters 
so much nor does our results at all. For the more accurate computations, in which we include all 
of the hadronic decay modes at one time, we would like to discuss them in separate issues. 

^Recently much milder value of the observational ^He mass fraction was reported by Olive and 
Skillman, 1^(05*) = 0.249± 0.009 [35]. Because their method is not the standard way to obtain 
the temperature and the electron density, we adopt it as a reference value here. The constraint 
which comes from Yp{OS) is denoted by the thin solid line in the figures hereafter. 
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gravitino decay becomes the dark matter of the universe. Fig. 3 depicts various constraints 
on the m^i2-m^ plane, (i) BBN constraints (thick and thin sohd hncs), (ii) Tr < mLSp/30 
(long dashed line), (iii) Tr > 1.2 MeV (dashed line), and (iv) > 2my2 (dot-dashed 
hne) in the case C = 10~^ and = 1. Also depicted there are two strips where neutrahno 
abundance from decaying gravitinos is in accord with the recently observed value Qlsp^^ = 
0.0945 - 0.1287 (95% C.L.) with mLSP = 100 GeV. The upper strip corresponds to the case 
gravitinos are thermally produced after modiihis decay. However, it is entirely covered by the 
region with Tr > mLSp/30 where the reheat temperature is so high that LSPs are thermally 
populated after modulus decay. In this case dark matter LSPs tend to be overpopulated 
unless some enhancement mechanism of annihilation is operative. This long dashed line has 
a scaling law for the other values of mLSP as oc mf^^gp. The lower strip, on the other hand, 
represents the case gravitinos are produced directly by modulus decay. There is a region on 
this strip where all the constraints are satisfied simultaneously with 777.3/2 — 55 — 100 TeV 
and 777<^ ~ (2 — 4) X 10^ TeV. As is seen in the figure it is highly nontrivial that such a region 
really exists. Since N appears only in the form N^/'^m^^/'^ in all of the relevant expressions 
in Eqs. (3), (10) and (14), the constraints for N other than N = 1 can easily be read off by 
replacing 777^ in the vertical axis by N^^^m^. Concerning this lower strip, there is a scaling 
relation 777,^ oc 777 lsp- Thus, for the other values of 777lsp we can easily see how the results 
change by using this relation. 

Fig. 4 depicts how the results change as we take different values of C which parameterizes 
the branching ratio of the modulus decay into the gravitino pair. We plot the cases (a) 
C = 10-^ (b) C = 10-^ (c) C = 10-2, and (d) C = 10-\ respectively From this panel, we 
can find that wide ranges of 777,^ are allowed (~ (9(10^) - O{10^) TeV) when C is changed for 
~ O(10~^) - C(10~^). Here, note that the intersection between the trip and the boundary 
of the BBN constraint gives a same value of 777,^ (~ 55 TeV) for any values of C. 

In Fig. 5 we depict the domain in the 7773/2-777^ plane where Qlsp^^ can take the observed 
value with some value of C in the range C ~ 10~^ — 10~^. The vertical thick (thin) solid 
line represents the BBN constraint, which comes from the observational '^He mass fraction 
by Izotov and Thuan Yp{IT) ( Olive and Skillman Yp{OS)). The long dashed line represents 
the boundary of Tr < 777lsp/30 . The dot-dashed line denotes the boundary of 777^ > 27773/2. 

As is seen there, we have a fairly large allowed region for 7773/2 ~ 55 - 2 x 10^ TeV and 
7770 ~ 10^ - 4 X 10^ TeV, with a small hierarchy between 777,^ and 7773/2. 

For different values of 7r7LSPi because the boundary of Tr < 7r7LSp/30 has a scaling law 

2 /3 

as oc 7r7Lgp, we can rescale the long dashed line by using this relation. In addition, about 
the modification of the vertical thick and thin lines for the other values of 777lsp, we can 
read off their values of 7773/2 as the intersections between the trip and the boundaries of the 
BBN constraint in Fig. 1 or Fig. 3. The long dashed line for N other than N — 1 can easily 
be read off by replacing 777,^ in the vertical axis by N^/^m^. 

VI. DISCUSSION 

We have shown that the neutralino dark mater whose abundance agrees with the WMAP 
constraint can be non-thermally produced by the decay of heavy gravitino with the mass 
of 7773/2 ~ 55 - 100 TeV in the early universe. Such gravitinos are produced only by the 



11 



decaying moduli with the late-time large entropy production, which completely dilutes both 
the thermal relics of LSPs, and primordial gravitinos produced by the reheating process 
after the primordial inflation. In this case, we found that the required region of the moduli 
mass is ~ 10^ - 4 x 10^ TeV. This scenario can solve the overproduction problem of 
LSPs in the supersymmetric cosmology. 

Note that although the LSP is still the best candidate for cold dark matter in particle 
physics, the recent measurement of the abundance by the WMAP restricts the preferred 
regions of the SUSY parameter space. The old and naive argument that the dark matter 
abundance is of the order of the critical density of the universe, which allows the wide 
regions of the parameter space, does no longer apply today because its abundance has been 
measured to be only one quarter of the critical density, which turned out to be significantly 
smaller than we expected in the standard calculation. As mentioned in the introduction, 
attempts do exist to explain this smallness within the standard thermal history, but this 
in turn requires us to adopt some enhancement mechanism of pair annihilation such as 
coannihilation or resonance enhancement etc (For the details see, e.g. recent works [36,37]). 
Our scenario is in a sense more economical than these mechanisms because we can account 
for the desired abundance of the neutralino dark matter as these particles are produced from 
modulus decay which also produces large entropy to solve the new gravitino problem at the 
same time. 

It is interesting to point out an alternative application of the production mechanism de- 
scribed in this paper. The neutralino production by gravitino decay can allow the neutralino 
dark matter with large fraction of the Wino and Higgsino components: its thermal abun- 
dance would be too small to constitute the dark matter of the universe. In this case, the 
cross section of the neutralino-nuclei scattering and the neutralino-neutralino pair annihila- 
tion can become larger. Thus, in the direct searches such as CDMS [38], EDELWEISS [39], 
and DAMA [40], we expect a higher detection rate than that of pure Bino LSP. In the 
indirect searches, we can also expect higher detectability of the signal of photons, neutri- 
nos, electrons and their radio emission through the neutralino-neutralino pair annihilation 
in the galactic center in the Milky Way, the center of the other galaxies, and the clusters of 
galaxies in the future experiments such as GLAST [41] and MAGIC [42] for gamma-rays, 
and Hyper-Kamiokande [43] and TITAND [44] for neutrinos, and so on. 

Finally we would like to emphasize that in this scenario the large gravitino mass will make 
the superparticle mass spectrum very different from the conventional gravity mediated SUSY 
breaking, and we expect the admixture of gravity mediation and anomaly mediation [45-48], 
which will be testable in future collider experiments. 
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FIGURES 




FIG. 1. Upper bound on the yield variable 7/3/2 = '^3/2/5- For heavy gravitinos with their 
masses of rn-3/2 > lOTeV, the constraint mainly comes from the ^He overproduction. The thick 
solid line, which is denoted by "IT", comes from the observational value of ^He mass fraction 
by Izotov and Thuan (2003). The dotted strip indicates the region where the LSP abundance 
from the gravitino decay corresponds to the desired value Olsp^^ = 0.0945 — 0.1287 (95% C.L.). 
The thin solid line, which is denoted by "OS", is the case that we adopted the most conservative 
observational value of Yp by Olive and Skillman (2004). Here we assumed that a gravitino can 
decay into all of the species of the chiral fermion and its scalar partner pairs (V'V') the gauge 
boson and the gaugino pairs (gg). 
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FIG. 2. Allowed regions for suitable observational value oin^sph^ = 0.0945-0.1287 (95% C.L.) 
in the C-m^ plane. We plotted the strips for m^i2 = 50, 10^, 5x 10^, and 10^ TeV, respectively. We 
considered only two conditions that Tr < mLSp/30 (below the long dashed line), and > 2m^/2- 
Here we took = 1, and rriLSP = 100 GeV. 
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FIG. 3. Various constraints on the m^/2^''^<f> plane, (i) BBN constraints (thick and thin 
sohd hnes), (ii) Tr < mLSp/30 (long dashed line), (iii) Tr > 1.2 MeV (dashed line), and 
(iv) > 2m3/2 (dot-dashed line) in the case C = 10~^ and = 1. Also depicted there are 
two strips where neutralino abundance from decaying gravitinos is in accord with the recently 
observed value Olsp/i^ = 0.0945 - 0.1287 (95% C.L.) with m-LSP = 100 GeV. There is a region 
on this strip where all the constraints are satisfied simultaneously with — 55 — 100 TeV 

and 7710 ~ (2 — 4) X 10^ TeV. Since appears only in the form N^l'^m^l'^ in all of the relevant 
expressions in Eqs. (3), (10) and (14), the constraints for N other than = 1 can easily be read 
off by replacing in the vertical axis by N^/^m^. 



17 



10 102 lono 102 




FIG. 4. How the results change as we take different values of C, which parameterizes the 
branching ratio of the modulus decay into the gravitino pair. We plot the cases (a) C = 10~^, 
(b) C = 10"'^, (c) C = 10~^, and (d) C = 10~^, respectively. For the explanation of the various 
constraints, see the figure caption in Fig. 3. 
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FIG. 5. Allowed region in the m^i2-Tn^ plane where JIlsp^^ can take the observed value with 
C ~ 10~^ — 10^^. The vertical thick (thin) solid line represents the BBN constraint, which 
comes from the observational ^He mass fraction by Izotov and Thuan Yp(IT) ( Olive and Skillman 
Yp{OS)). The long dashed line represents the boundary of Tr < ?tilsp/30 . This line for N other 
than = 1 can easily be read off by replacing in the vertical axis by N^/^m^. The dot-dashed 
line denotes the boundary of m,^ > 2777.3/2 • 
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